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DECLARATION UNDER 37 C.F.R, § 1.131 



Commissioner for Patents 

United States Patent and Trademark Office 

P.O. Box 1450 

Alexandria, VA 22313-1450 

Dear Sir 

1. This declaration is to establish the status of the invention in the above- 
captioned U.S. Patent Application prior to March 26, 1993, the filing date of U.S. Patent 
No. 5,368,026, issued on November 29, 1994, to Swedlow et al. ("Swedlow"). 

2. We understand that Claims 1-10, 12-16, 18-25, 27-36, and 38-46 of the 
above-captioned U.S. Patent Application Ser. No. 10/677,050 (the '050 Patent 
Application), which was filed September 30, 2003, titled "Signal Processing Apparatus," 
have been rejected in the February 27, 2008 Office Action as being unpatentable over 
Swedlow. 

3. We are inventors of the subject matter of the above-referenced claims in 
the '050 Patent Application. 

4. At least prior to March 26, 1993, we conceived of and reduced to practice 
a signal processing apparatus such as the device and method described and claimed in 
the above-referenced claims of the '050 patent application, 

5. As evidence thereof, we provide herewith Exhibit A which is a design 
note ("the Design Note") describing the subject matter of the above-referenced claims. 
The Design Note was created and is dated prior to March 26, 1993. We conceived the 
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ideas referenced therein prior to the creation of the Design Note. Any references to 
dates on the Design Note have been redacted. 

6. At Masimo Corporation ("Masimo"), the term "Stork" refers generally to 
Masimo's pulse oximetry technologies. 

7. At least prior to March 26, 1993, we compiled, loaded onto hardware, and 
used a software application called "STORK" for its intended purpose, that is, to measure 
blood oxygen saturation. 

8. As evidence thereof, we provide herewith Exhibits B and C which are 
screen captures of a directory of files, including "\c" source code files, ".h" header code 
files, and ".gdd" debugging code files compiled into the application file STORK, As 
shown on the screen captures, each of these files includes a last edited date prior to 
March 26, 1993, which has been redacted herein. 

9. Exhibits D and E are the header and source code files, titled "vector.h," 
and Vecfor.c," respectively, written in the C programming language, showing the 
mathematical functions for calculating a vector cross correlation. 

10. Exhibits F and G are the header and debugging files, titled "stork.h" and 
"stork.gdd," respectively, written in the C programming language. The header file 
defines the variable for cross correlation and the debugging file displays the value of 
cross correlation for instructive purposes. 

11. Exhibit H is the source code file, titled "iawa.c," written in the C 
programming language, which calculated the cross correlations. 

12. Exhibit I is the source code file, titled "satfil.c," written in the C 
programming language, which qualifies bins of signal data responsive to intensity 
signals from a detector monitoring a patient, the qualification including determining 
whether the calculated cross correlation is within expected limits. 

13. Exhibit J is the source code file, titled "ui.c" written in the C programming 
language, which is the user interface for the STORK application and which indicates an 
expected limit of cross correlation. 

14. The foregoing files and other exhibits demonstrate that we constructed 
and tested an embodiment and performed processes meeting every element of the 
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above-referenced claims prior to March 26, 1993 sufficient to demonstrate that the 
constructed embodiment operated for its intended purpose. 

15. We hereby declare that acts leading to the reduction to practice of the 
subject matter claimed in the '050 Patent Application were performed in the United 
States. 

I declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true, and further that 
these statements were made with the knowledge that willful, false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful, false statements may jeopardize 
the validity of the application or any patent resulting therefrom. 



I declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true, and further that 
these statements were made with the knowledge that willful, false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful, false statements may jeopardize 
the validity of the application or any patent resulting therefrom. 



Dated; By: 



Dated: TTua^c 13 - t 




Mohamed K. Diab 



Walter M. Weber 
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I declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true, and further that 
these statements were made with the knowledge that willful, false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful, false statements may jeopardize 
the validity of the application or any patent resulting therefrom. 



By: 

Massi E. Kiani 
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If 


WIIKin CORRELATION DETECTION n^ 9 " 
MU ™ ALGORITHM n ° iES 


Design Nolc 
37 


Revision 

A 


File 

SDN37A.DOC 


Stork 


Page 

lof2 


Walter Weber 


Revision Dale 


Original Dale 







The goal of this note is to propose an artifact detector for the Stork algorithm. Specifically, the artifact 
detector is designed such that it assigns a real number to each snapshot of data that it processes, and then 
compares that real number to a predetermined threshold value. Snapshots of data which yield artifact 
detector values that lie below the threshold value are said to contain significant amounts of artifact. 
Snapshots of data that yield artifact detector values which lie above the threshold value are said to contain 
insignificant amounts of artifact. 



In principle, the Stork algorithm accepts as input two signals which a 
bandpass) 



s a(0 = s a(0 + »a(0 
S b (0 - s D («) + n bW 



; of the form (after log and 



Any input signals which are not ex 
Consequently to determine whether or i 
not equations (1) and (2) are satisfied. 



s b(0 = »1 s a(0 



ssable in the above form are said to consist of artifacts, 
a signal is artifact free, it is sufficient to determine whether or 



-"One way to check whether or not equations (1) and (2) have been met is to n 
between S a (f) and S b (l). Specifically one could compute the cross correlation 



f S a (t) S b (t) dt = I ( s a (t) s b (t) + n a (t) n b (t) ) dt + I ( s a (t) n b (t) + s^t) n a (t) ) dt 

that the signals, s a (t) and s b (t), are uncorrected with the noises, n a (t) and n b (l), may be 



J S a (0 S b (t) dt = J( s a (t) s b (t) + n a (t) n b (t) ) dt. 

A normalized cross correlation may be computed as follows 

J S a (t) S b (t) dt I ( s a (t) sb<t) + n^t) n b (t) ) dt 



fj S 2 a (t) dt I S 2 b (t) dt]°-5 [J ( s 2 a(t) + „2 a(t) ) dt ( ( s 2 b(t) + n 2 b(t) ) dtJ 0.5 



The solutions which give the extrema of this cross correlation functional are given by equation (2). Any 
input signals which are not expressible in the form given in equation (2) yield reduced normalized 
correlation values. In Figure 1, these extrema of the normalized correlation functional are plotted vs. the 

Copyright ©^■Masimo Corporation ™ 26052 Merit Circle Suite 103, Laguna Hills, CA, 92653 
^^This document contains confidential and proprietary information. 



"Red" Signal - to - Noise ratio and m x lm 2 . One should note that the normalized correlation is not 
necessarily unity and precaution must be~taken, particularly when the difference between arterial 
saturation and the background saturation is large. 

One may rigorously prove that the magnitude of the normalized cross correlation is always less than or 
equal to unity by use of the inner product 

x»J =fx(t)y<t)dt. 
Using the above notation, the normalized cross correlation may be written as 

V«b + »a»% 
' KVSa + n a*"a)(Sb*s b + %»»b)]°- 5 

or equivalently 

Sa'Sp 

IVSaVSb] 0 - 5 ' 

But, according to the Schwarz Inequality 

(S a .S b )2 < (S a .S a )(S b .S b ). 
Hence the normalized correlation is always between 1 and -1. 
ALGORITHM IMPLEMENTATION 

1 . Compute the quantity 

£s a s b 

IzPTz s2 b] <>-5 " 

2. Compare quantity to a threshold, say 0.75 

THINGS THAT STILL NEED TO BE INVESTIGATED 

Investigate performance vs. snapshot size 
Investigate performance vs signal/noise correlation 
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Other Places 



©ADC.CC 
©HW.H 



3KB CHeaderffe 

4KB C Source Re 

4KB CHeaderrte 



6 KB C Source fie 

2 KB c Header fie 

3 KB CSourcefte 
3KB CHeaderrte 
10 KB C SOITCB fie 

3 KB CHtederfte 

5 KB csourcefte 

4KB CHesderfte 

12KB CSourcefte 

7 KB CSourcefte 
2KB CHeadEffte 
4KB C source fte 
2KB CHeaderfte 
4KB CSourcefte 
4KB CSourcefte 
2KB CHearJa-fte 

13KB CSourcefte 



1KB CHeadH-rte 
4KB CHeaderfle 
2KB C Header fit 
4KB CSourcefte 
4KB CSourcefte 
7 KB Fife 
3 KB csourcefle 
6 KB CSourcefle 
3KB CHeaderrte 

9 KB GOD Re 

10KB CHeaderrte 

WKB CSourcefte 

10 KB C Source fit 

0 KB OAT Fife 

0 KB Text Document 

6KB CSourcefte 

30KB CSourcefte 

WKB CSourcefle 

260 KB A 



2A9A993 11:27 AM 
2/19/1993 1 1:40 AM 
2/13/1393 11:45AM 



2/23/1993 1:26 PM 
2/23/19933:15 PM 
2/23/1993 5:«PM 
2/23/1993 S:59PM 
2/24/1993 3:02 PM 
2/24/1993 3:03 PM 
2/25/1993 12:00 PM 
2/25/19931:00 PM 
2/2S/1993 4:26 PM 
2/25/1993 4:26 PM 
2/26/1993 5:S0PM 
3/2/1993 10:09 AM 
3/2/1993 l:12PH 
3/2/1993 2:45 PH 
3/2/1993 2:47 PH 
3/2/1993 2:47 PU 
3/2/1993 5:51PH 
3/3/1993 S:03PH 
3/3/1993 5:03 PH 
3/3/1993 S:2BPH 
3/S/lB93 2:2lPM 

3/5/1993 2:54 PH 
3/5/1993 3:06 PH 
3/5/1993 3:06 PH 
3/5/1993 3:27 PH 
3/5/1993 3:37 PH 
3/5/1993 1:36 PH 
3/6/1993 IWIPH 
3/8/1993 8:58 AM 
3/8/1993 10:21 AM 
3/8/1993 12:47 PM 
3/8/1993 12:50 PM 
3/8/1993 1:04 PH 
3/B/1S93 1:07 PH 
3/8/1993 1:15 PH 
3/8/1993 2:06 PH 
3/8/3993 3:44PH 
3/6/1993 3:45 PH 
3/B/1993 3:4B PH 
3/8/1993 3I49 PH 
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©ASH.SPRT.H 
@0£F21020.H 
©HANOI 

©VSIMCM 
©STATS.H 
©PEAK.H 
©SATFILH 
©HRHLH 



©NRLCOfUl 

©HaiXOR-C 

®H(LH 

©ZFRJ^C 

©F05TOS_CC 

®ZPOSTC6_.C 

©OPTANCH 

©CDF.H 

©ZUSftl_CC 

©CDF.C 



3KB CHtadcrrte 

8KB CBeaderfte 

1KB CHeaderfle 

1KB c Header He 

2KB CHeadH-fte 

1KB C Header rit 

2KB CHeadErfte 

1KB CHtftdtrflt 

2KB CHeada-fte 

2KB CHeaderfte 

2KB C Header fie 
ndentramarifr derffe 

7 KB CSourcefle 



3 KB CSourcefb 

3KB CHeaderfte 

3 KB CSourcefle 

8KB CSourcefle 
8KB CSourcefle 
1KB CHeaderfte 
2KB CHeaderfte 

6 KB CSourcefle 

4KB CSourcefte 



Tcwe Mod fad + 



©Z1CMJBX 

©zivsmcc 

©Z1HANC.C 

©Z1SA7HLC 

©Z1STATS.C 

©ZIOPTANCC 

©ZiVSMRjC 

©Z1VSIPLC 

©Z1HRHLC 

©Z1C0F.C 

©ZllSLC 

©Z1SATKJN.C 

©Z1IA.WA.C 

©Z1A0CC 

©Z1PACKET.C 

@ZlHa.COR.C 



1KB CSourcefle 

2KB CSourcefte 

3 KB CSourcefle 

1KB CSourcefle 

1 KB CSourcefle 

1KB CSourcefle 

1KB CSourcefle 

1KB CSourcefle 

2KB CSourcefle 

1KB C Source Re 

1KB csourcefle 

1KB CSourcefte 



If 12/1932 2:43 PM 
7/21/1992 11:27 AM 
11/6/1992 5:33 PM 
11/10/1992 11:29 m 
12/8/1992 3:55 PM 
12/30/1992 4:53 PH 
12/30/1992 4:34 PM 
l/t6/1993 2:34 PM 
1/18/1993 2:38 PM 
1/18/19935:11 PM 
1/25/1993 2:14 PM 
2/3/1993 4:29 PM 
3/4/1993 11:22 AM 
2/9/1993 2:23 PM 
2/9/1993 2:24 PM 
2/9/1993 2:26 PM 
2/12/19934:40 PM 
2/12/19934:52 PM 
2/16/1993237 PM 
2/16/1993 2 ;38PM 
2/16A9932:4SfM 

2/16/1993 2:45 PM 
2/t6/19932:50PM 
2/17/1993 10:32AM 



2/17/1993 2 :30 PM 
2/18/19933:11 PM 
2/lB/19335:45PM 
2/W19936:S3PM 
2/18/1993 6:56 PM 
2/18/19936:57 PM 
2/18/19936:58 PM 
2/18/1993 6:59 PM 
2/18/1993 7 :02 PM 
2/18/19937:04 PM 
2/18/1993 7 :06 PM 
2/18/1993 7 :08 PM 
2/18/19937:10 PM 
2/18/19937113 PM 
2/18/19937:17 PM 
2/16/19937:19 PM 
2/18/1993 7 :20 PM 
2/lB/1993 7:22PM 
2/18/1993 7:26 PM 
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1 /* 

Copyright (c) Masimo Corporation (tm) 1993. All Rights Reserved. 



File: vector. h 

Description: 

Collection of routines that work on vectors. 
Masimo Standard Library. 

Public Functions: 
min_index 
max_index 
inner_jproduc t 
remove_mean 



Returns index to minimum element in vector 
Returns index to maximum element in vector 
Returns point by point product of vectors 
Removes the mean value from each point in a vector 



VECTOR_DotProduct 
VECTOR_CrossCorrO 



Private Functions : 



Computes Dot product of two vectors 
Computes Correlation Coefficient 



01/15/93 
01/25/93 
02/25/93 
02/25/93 



RJM 
RJM 
RJM 
RJM 



Initial Version 
add remove_mean 
add VECTOR_DotProduct 
Add VECTOR_Cros sCore 0 



#include <masimo.h> 



extern INT32 
min_index ( 
FLOAT32 
INT32 



*data, 
size) ; 



extern INT32 

max_index ( 

FLOAT32 *data, 

INT32 size) ; 

extern VOID 

innerjproduct ( 

FLOAT32 *vl, 

FLOAT32 *v2, 

INT32 size, 

FLOAT32 *v3); 

extern VOID 

remove_mean ( 

FLOAT32 *in, 

FLOAT 3 2 *out, 

INT32 size) ; 

extern FLOAT32 
VECTOR_Dot Product ( 

FLOAT32 *u, 

FLOAT 3 2 *v, 

INT32 size) ; 

extern FLOAT 3 2 
VECTOR_CrossCorrO ( 

FLOAT 3 2 *u, 

FLOAT 3 2 *v, 

INT32 size) ; 



/* input, pointer to data 

/* input, number of data points 



input, pointer to data 
input, number of data points 



/* input, first vector 

/* input, second vector 

/* input, size of vectors 

/* output, v3 = vl * v2 



/* input, data 

/* output, data with mean removed 
/* input, number of points 



/* input, first vector 
/* input, second vector 
/* input, size of vectors 



/* input, first vector 
/ * input , second vector 
/* input, size of vectors u and ■ 



Exhibit E 



1 /* 

Copyright (c) Masimo Corporation (tm) 1993. All Rights Reserved. 



File: vector. c 

Description: 

Collection of routines that work on vectors. 
Masimo Standard Library. 



Public Functions: 
min_index 
max_index 
innerjproduct 
remove_mean 
VECTOR_DotProduct 
VECTOR_CrossCorrO 



Designed for 



Returns index to smallest element in vector 
Returns index to largest element in vector 
Computes the inner product of two vectors 
Removes the mean value from each point in a i 
Compute Dot Product of two vectors 
computes cross correlation time=0 



Private Functions: 

Notes: 

History: 

01/15/93 RJM Initial Version 

02/25/93 RJM Add VECTOR_DotProduct 

02/25/93 RJM Add VECTOR_CrossCorrO 



#define MODULE_ID 

#include <masimo.h> 
#include <math.h> 
#include "vector. h" 



/* platform definitions 
/* sqrt function 
/* self 



size)"' { 



INT32 

min_index ( 
FLOAT32 
INT32 

INT32 
INT32 
FLOAT 3 2 



BUG1 (data) ; BUG1 (size) ; 



/* input, pointer to data 

/* input, number of data points 



tmp = data [ j ] ; 
for (i = 1; i < size; i++) 
if (data[i] < tmp) { 
tmp = data [i] ; 



INT32 

max_index ( 

FLOAT32 *data, /* input, pointer to data */ 

INT32 size) { /* input, number of data points */ 

INT32 i; 

INT32 j ; 

FLOAT32 tmp; 



« 

BUG1 (data) ; BUG1 (size) ; 



j = 0; 

tmp = data [ j ] ; 

for (i = l; i < size; i++) { 
if (data[i] > tmp) { 
75 tmp = data [i] ,- 



return (j ) ; 

80 



VOID 

inner_product ( 
FLOAT32 
FLOAT32 
INT32 
FLOAT 3 2 

INT32 



for (i = ( 
V3 [i] 

} 



*vl, 
*v2, 
size, 
*v3) { 

i; 



/* input, first vector 

/* input, second vector 

/* input, size of vectors 

/* output, v3 = vl * v2 



i < size; i++) { 
vl [i] * v2 [i] ; 



VOID 

remove_mean ( 
FLOAT 3 2 
FLOAT32. 
INT32 



*out, /■ 
size) { /' 



input, data 

output, data with mean removed 
input, number of points 



BUG1 (in) ; BUG1 (out) ; BUG1 (size) ; 
mean = 0.0; 

for (i = 0; i < size; i++) mean += in[i] ; 
mean /= (FLOAT32) size; 

for (i = 0; i < size; i++) out [i] = in[i] 



/* sum all points */ 
mean; 



FLOAT 3 2 

VECTOR_DotProduct ( 
120 FLOAT32 *u, /* 

FLOAT32 *v, /* 

INT32 size) { /* 

FLOAT32 sum; 
125 INT32 i; 

BUG1 (u) ; BUG1 (v) ; BUG1 (size) ; 

sum = 0.0; 
130 for (i = 0; i < size; i++) { 

sum += u[i] * v[i] ; 

} 

return (sum) ; 



input, first vector */ 
input, second vector */ 
input, size of vectors */ 



135 



} 



FLOAT32 

VECTOR_CrossCorrO ( 
FLOAT32 
FLOAT32 
INT32 



*u, 



/* input, first vector 
/* input, second vector 
/* input, size of vectors u and v 



BUG1 (u) ; BUG1 (v) ; BUG1 (size) ; 

den = SQRTF (VECTOR_DotProduct (u, u, size) * 
VECTOR_DotProduct (v, v, size) ) ; 

if (den ==0.0) { 

c = 1.0; 
} else { 

c = VECTOR_DotProduct (u, v, size) / den; 

} 

return (c) ; 
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Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 
File: stork. h 

Description: Stork header file included by ALL stork modules 
Public Functions: 



10 Private Functions: 



Notes: 



History: 




RJM 


11/06/92 


RJM 


12/08/92 


RJM 


12/16/92 


RJM 


01/04/93 


RJM 


01/05/93 


RJM 


01/14/93 


RJM 


01/15/93 


RJM 


01/16/93 


RJM 


01/18/93 


RJM 


01/18/93 


RJM 


01/19/93 


RJM 


01/25/93 


RJM 


01/25/93 



Initial version {vl.29) 
Prevent multiple inclusion of stork. h 
Change vsinc params to all floats 
Add MAX_BINS 

Add seed to bin structure 

Put conditional code around "C" only stuff 
Make memory pool sizes platform dependant 
Change sat result structure for filtering 
Make HR a vector 
Add more sat params 

Add hr_max_delta and sat_max_delta params 
Make SAT equation adjustable 
Make HR algorithms switchable 



1.541 



#define VERSION 
#define VSINC_LEAN 0 
#define ADC_SNAPSHOT SIZE 



/* floating point software version number */ 
/* Generate code for Lean Architecture */ 



Memory Allocation 



#if (PLATFORM = 



#define MEM_STACK_SIZE 48 

#define MEM_POOL_SIZE 32768 

ttdefine DEBUG_MEM_SIZE 512 

50 #define GDD_BUF_SIZE 2 000 0 

#define DBG_FIFO_SIZE 100 



#else 

55 #define MEM_STACK_SIZE 48 

#define MEM_POOL_SIZE 26000 

#define DEBUG_MEM_SIZE 512 

#define GDD_BUF_SIZE 6000 

#define DBG_FIFO_SIZE 100 

60 

#endif 



/* data memory */ 

/* data memory */ 

/* data memory */ 

/* pmda (was 6000) */ 

/* pmda */ 



/* The following contant controlls ALL debug code for the STORK */ 
/* beta software. No required code should be conditioned on this */ 
65 /* constant. */ 



# define DEBUG 1 



/* enables debug code 



*/ 



/' 



The following constants control what data is captured in the dbg_fifo. 
The dbg_fifo is controlled by the ADC interrupt routine. 



# define DEBUG_SIM 
#define DEBUG_ADC 
ttdefine DEBUG_LED 
#define DEBUG_SUMS 
ttdefine DEBUG_RAW 
#define DEBUG_HSL 



simulate ADC data */ 

stores raw adc conversions in fifo */ 

stores tick mark for LED change */ 

stores the reading for each state */ 

stores raw ir and raw red */ 

stores High Speed Link data */ 



85 Low Speed Link Control 

*/ 

#define LSL_CHECK_OVERRUN 1 /* check HC11 FIFO */ 

#define LSL_GDD_GROUP_SIZE 6 /* number of gdd packets per lsl_tick */ 

90 

/* 

High Speed Link Control 
. */ 

95 #define HSL_RECORD 0 /* records ADC data after ambient subtract */ 

ttdefine HSL_PLAYBACK 0 /* uses data from HSL instead of ADC */ 

#define HSL_TEST_PATTERN 0 /* generate rotating l's pattern */ 

#define IMPULSE 0 /* generate impulses for usrl */ 

100 

/* - - 

The following constants control the ANC. 

*/ 

105 

#define ANC_FIXED_CELL 0 /* fixed cell normalized anc */ 

/* 

110 System Wide Type definitions 



ttdefine FIXED_CELL 1.0 /* anc_types */ 

115 #define ACCORD IAN 2.0 /* anc_types */ 

/* The following constants describe the size of the following */ 

/* structures. */ 



120 #define MAX_BINS 6 

#define NUM_HR_ALG 7 

#define BIN_RESULTS_SIZE 8 

#define SAT_PARAMS_SIZE 22 

125 #define SAT_RESULTS_SIZE 14 

ttdefine HR_PARAMS_SIZE 7 

ttdefine HR_RESULTS_SIZE 8 

ttifndef ASSEMBLY 

130 

typedef struct { 



INT32 


num_jpoints; 


/* 


number of input points 


*/ 


FLOAT32 


anc_type; 


/* 


type of anc 


*/ 


FLOAT32 


lambda ,- 


/* 


nANC parameter 


*/ 


FLOAT32 


min_error ; 


/* 


nANC parameter 


*/ 


INT32 


anc_settle_time ; 


/* 


nANC parameter 


*/ 





FLOAT 3 2 


num_cells ; 


/* 


nANC parameter 


*/ 




FLOAT 3 2 


peak_limitl; 


/* 


peak picker parameter 


*/ 




FLOAT 3 2 


peak_epsilon; 


/* 


peak picker parameters 


*/ 


140 


FLOAT32 


*bp_coef ; 


/* 


pointer to filter coef . 


*/ 




INT32 


bp_size; 


/* 


number of filter coef . 


*/ 




FLOAT3 2 


search_f irst ; 


/* 


first saturation value. 


*/ 




FLOAT3 2 


search step; 


1* 


saturation step size. 


*/ 




INT32 


search_size; 


1* 


number of steps 


*/ 


145 


INT32 


bin_count ; 


/* 


number of bins 


*/ 




FLOAT3 2 


coyJLimit; 


/* 


covariance limit 


*/ 




FLOAT32 


DVHS_init ; 


/* 


initial value for DVHS 


*/ 




FLOAT 3 2 


DVHS0_init ; 


/* 


initial value for first DVHS 


*/ 




FLOAT3 2 


DVAST_timer ; 


1* 


seconds 


*/ 


150 


FLOAT 3 2 


DVAST_delta ; 


/* 


change for DVAST 


*/ 




FLOAT3 2 


DWST timer; 


/* 


seconds 


*/ 




FLOAT32 


DWST_delta; 


/* 


change in DWST 


*/ 




FLOAT 3 2 


max_sat_delta ; 


1* 


max change in sat w/motion 


*/ 




INT32 


sat ecpi; 


1* 


saturation equation number 


*/ 


155 


FLOAT3 2 


dual_lambda ; 


/* 


dual ANCs as VSINC1 inputs 


*/ 




FLOAT3 2 


correlation; 


/* 


bins OK only if >= corr. 


*/ 




} SAT_PARAMS ; 












typedef struct 


{ 








160 


FLOAT32 


ir_rms; 


/* 


signal rms 


*/ 




FLOAT 3 2 


rd_rms ; 


/* 


signal rms 


*/ 




FLOAT 3 2 




/* 


covariance 


*/ 




FLOAT32 


arterial sat ; 


/* 


resulting arterial sat 


*/ 




FLOAT32 


venous_sat ; 


/* 


resulting venous sat 


*/ 


165 


FLOAT32 


seed_sat ; 


/* 


seed saturation 


*/ 




FLOAT 3 2 


num_peaks ; 


/* 


peaks in power curve 


*/ 




FLOAT32 


correlation; 


/* 


raw dataset correlation 


*/ 




} BIN RESULTS; 











170 typedef struct { 



FLOAT 3 2 


Spa02_f inal; 


/* 


resulting arterial sat 


*/ 


FLOAT32 


Spv02_f inal; 


/* 


resulting venous sat 


*/ 


FLOAT 3 2 


Spa02 unfil; 


/* 


unfiltered arterial sat 


*/ 


FLOAT32 


Spa02_fill; 


/* 


first stage of filtering 


*/ 


FLOAT32 


Spa02 dual; 


/* 


dual ANCs as input to VSINC1 


*/ 


FLOAT32 


DVHS; 


/* 


default value high sat 


*/ 


FLOAT32 


DVHSO; 


/* 


default value high sat first 


*/ 


FLOAT32 


DVAST; 


/* 


default value art sat (PR) 


*/ 


FLOAT32 


DVAST time; 


/* 


time DVAST was last updated 


*/ 


FLOAT32 


DWST; 


/* 


default value ven sat (PR) 


*/ 


FLOAT32 


DWST_time; 


/* 


time DWST was last updated 


*/ 


FLOAT 3 2 


art_omega_pr ,- 


/* 


arterial omega for pulse rate 


*/ 


FLOAT32 


ven_omega_pr ; 


/* 


venous omega for pulse rate 


*/ 


FLOAT32 




/* 


time for these results 


*/ 


BIN RESULTS 


bins [MAX_BINS] ,- 


/* 


bin results 


*/ 



} SAT_RESULTS ; 



typedef struct { 








FLOAT 3 2 


resolution; 


/* 


HR resolution in BPM 


*/ 


FLOAT 3 2 


raw_sampling_rate ; 


/* 


raw data sampling rate in Hz 


*/ 


FLOAT 3 2 


min hr; 


/* 


minimum hr allowed (BPM) 


*/ 


FLOAT32 


max_hr; 


/* 


maximum hr allowed (BPM) 


*/ 


FLOAT32 


nst; 


/* 


noise spec thresshold 


*/ 


FLOAT32 


max_hr_delta ; 


/* 


maximum change in HR w/motion*/ 


INT32 


hr_enable; 


/* 


bits control with alg 


*/ 


FLOAT 3 2 


hr_lambda ; 


/* 


ANC algo. uses this lambda 


*/ 


} HR PARAMS; 










typedef struct { 








FLOAT32 


hr [NUM_HR_ALG] ; 


/* 


(BPM) from each alg 


*/ 


FLOAT32 


hr_f inal; 


/* 


(BPM) after all filtering 


*/ 


FLOAT32 


BPM_jper_bin; 


/* 


BPM per bin of spectrum 


*/ 


INT32 


num_bins ; 


/* 


number of bins per spectrum 


*/ 



205 ^} HRJRESULTS; 
#endif 
#endif 



Exhibit G 



Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 
File: stork. gdd 

Description: Debug Display Control File 
Public Functions: 
Private Functions: 
Notes : 

Channels 0 to 39 are defined in SPT.H 



15 


Channels 97 to 88 


are for status line j 




History: 








RJM 


12/10/92 


Initial version 


20 


RJM 


12/21/92 


Add algorithm parameter \ 




RJM 


12/22/92 


Add bin results table and plots f 




RJM 


12/28/92 


Add new HR algorithm results 




RJM 


12/30/92 


Change to one algorithm J 




RJM 


12/30/92 


Add peak_epsilon parmaters 


25 


RJM 


01/05/93 


Add Seed sat to bin data, remove others 




RJM 


01/16/93 


New SAT results structure ] 




RJM 


01/18/93 


Add dvhs_init and dvhs0_init parameters 




RJM 


01/18/93 


Add hr5 result 




RJM 


01/18/93 


Add more sat params 


30 


RJM 


02/22/93 


Real-Time data Plotting 




RJM 


02/23/93 


Clean up Channel Assignments 


35 


Debugl 


screen 





P 0,2 
DEBUG 
30,050, 
34, 050, 
33,050, 
32,050, 
20,775, 
21,775, 
22,775, 
23,775, 
24,775, 
25, 775, 
31,775, 
00,400, 
01,575, 
03,575, 
04,575, 
02, 010, 
17,400, 
98, 000, 
99,000, 
97,000, 



025,0375 
175, 0375. 
325, 0375 
475, 0375 
025, 1000 
150,1000 
230,1000. 
310, 1000, 
390, 1000, 
470,1000, 
600, 1000, 
005, 0000, 
005, 0000, 
225, 0000, 
370, 0000, 
610,0000, 
490, 0014, 
741,0014 
741, 0015, 
755, 0015, 



125,7 
275,5 
425,5 
575,5 
100,5 
225,0 
305,0 
385, 0 
465,0 
545, 0 
725,5 
000, 0 
000,0 
000, 0 
000,0 
000, 0 
000,12 
000, 12 
000,12 
000, 12 



; display key = : 

; raw data 

; analysis data 

; plethy 

; freq spect 

; power curve 

; bin 1-5 



ADC Snapshot 
system params 
Sat Results 
HR Results 
Calib Results 
bin results 
ADC OOR count 
' Working ' Text 
Version# 
file name . . 



; Debug2 screen 

V 1,1 

DEBUG SCREEN 2 
30,050,025,425,150,5 
32, 050, 225,425, 330, 7 
33,050,400,425,525,5 
00,450, 005, 000, 000, 0 
01,625,005,000,000,0 



Raw, Spectrum, Plethy, tables 
; display key = 2 
; raw 
; plethy 

; System Parameters 
,- Sat Results 



03,625,225,0000,000,0 
04,625,400,0000,000, 0 
98,050,570,014,000,12 
99,050,570,015,000,12 
97,050,585,015,000,12 



HR Results 
Calib Results 
'Working' Text 
Version# 
file name. . 



Debug3 screen power curve, bin graphs 



P 2,1 

DEBUG SCREEN 3 

20,050,020,380,130,7 

21,050,180,380,280,7 

22,050,330,380,430,7 

23,450,020,780,130,7 

24,450,180,780,280,7 

25,450,330,780,430,7 

02,060,460,000,000,0 

98,050,590,014,000,12 

99,050,590,015,000,12 

97,350,590,015,000,12 

END 



; display key = 3 



bin results 
1 Working ' Text 
Version# 
file name. . 



Debug4 screen Raw, Spectrum, Plethy, tables 



P 3,0 
DEBUG ; 
30,050 
32,050 
33,050 
00,300 
01,470 
03,470 
04,470 
98, 000 
99,000 
97, 000 



SCREEN 4 
025,270,125,5 
180, 270,280,5 
330,270,430,4 
,-15,000,000,0 
-15,000,000,0 
195, 0000, 000, 0 
320,0000,000,0 
,455,014,000,12 
,455,015,000,12 
,470,015,000,12 



display key = 4 



plethy 

system params 
sat results 
HR Results 
Calib Results 
'Working' Text 
Version# 
file name. . 



Debug5 screen power curve, bin graphs 



P 4,0 

DEBUG SCREEN 5 

20,050,020,295,100,7 

21,050,135,295,215,7 

22,050,250,295,320,7 

23,375,020,620,100,7 

24,375,135, 620,215,7 

25,375,250, 620,320,7 

02,020,345,000,000,0 

98,020,471,014,000,12 

99,020,471,015,000,12 

97,250,471,015,000,12 



bin results 
' Working 1 Text 
Version# 
file name. . 



; Customer Screen #1 

P 5,0 

Customer Screen 

52,000,000,015,000,28 

77,050,050,450,150,21 

78,050,200,450,300,21 

53,500,080,014,000,16 

54,500,230,014,000,16 

15,500,110,014,000,16 

16,500,260,014,000,16 

31,500,310, 630,400,00 

17,500,415, 011, 000, 08 



display key = 6 

stork legend 

Real time IR Data (23) 

Real time RD Data 

Spa02 : 

Pulse : 

Spa02 value 

Pulse value 

ADC Snapshot 

ADC OOR count HGK 



18,500,430,011,000,0 
14,500,445,011,000, 0 
40,500,460,012,000,0 



Correlation 
Tick 

calibrate: 1/0 message 



Customer Screen #2 

P 6,0 

Customer Screen 

52,000,000,015,000,28 

35,050,050,450,150,21 

33,050,200,450,300,05 

53,500,080,014,000,16 

54, 500,23 0,014, 000, 16 

15,500,110,014,000,16 

16, 500,260,014, 000, 16 

31, 500,310, 630,400, 00 

17,500,415,011,000,08 

18,500,430,011,000,08 

14, 500,445, 011, 000, 08 

40,500,460,012,000,08 



display key = 7 

stork legend 

raw analysis portion 

plethy 

Spa02 : 

Spa02 value 

Pulse value 

ADC Snapshot 

ADC OOR count HGK 

Correlation 

Tick 

calibrate: 1/0 message 



; Front End Screen 

P 7,1 

Front End Screen 

04, 500, 010, 0000, 000, 0 

31,500,200,775,575,7 

17,050,050,014,000,12 

77,050,200,450,375,23 

78,050,400,450,575,23 

40,050,020,011,000,28 



Calib Results 
ADC Snapshot 
ADC OOR count 
Real time IR Data 
Real time RD Data 
calibrate: 1/0 message 



Channel Definitions 



Usrl snapshot 



G 30 

Raw Data 
Point number 
Amplitude 



ADC snapshot 



G 31 

ADC SNAPSHOT 
Point Number 
Amplitude 



Frequency Spectrum 



G 32 

Spectrum 
HR (BPM) 
Amplitude 



Master Power Curve 



G 20 

Power Curve 
Sat 



G 21 
Bin #1 
Sat 



G 22 
Bin #2 
Sat 
Amp 



G 23 
Bin #3 
Sat 
Amp 



G 24 
Bin #4 
Sat 
Amp 



G 25 
Bin #5 
Sat 
Amp 

Plethysmograph 

G 33 

Clean Plethysmograph 

Time 

Amp 

; analysis data 

G 34 

Analysis Data 
Amp 

; raw analysis data 

G 35 

Raw Data (before Masimo) 

Index 

Amp 

; Info Strings 

F 98,0 

Working 

F 99,0 

Stork Version %6.3f 

; Algorithm Parameters 



T 0, 1, 32 
ALG 1 

covJLimit, %9.1f 
num_jpoints , %9 . Of 



anc_type, %9.0f 
lambda, %9.4f 
min_error, %9 .4f 
anc_settle,%9.0f 
num_cells,%9.0f 
peak_limitl, %9.2f 
peak_epsilon, %9 . 2f 
search_f irst , %9 . 2f 
search_step, %9 . If 
search_size, %9 . Of 
bin_count,%9.0f 
sampling rate, %9 . Of 
dvhs_init,%9.1f 
dvhsO_init , %9 . If 
DVAST_timer, %9 . Of 
DVAST_delta,%9.1f 
DWST_timer , %9 . Of 
DWST_delta, %9 . If 
max_sat_delta , %9 . If 
max_hr_delta , %9 . If 
sat eqn,%9.0f 
hr enable, %9. Of 
Cal Target, % 9. 4 f 
Cal Err Band, %9.4f 
Cal Set. Time, %9. Of 
Cal IR PGA, %9.2f 
Cal RD PGA, %9.2f 
hr lambda, %9.4f 
dual lambda , % 9 . 4 f 
Correlation, %9.4f 

Sat Results 

T 1,1,14 
ALG 1 

>Spa02_f inal , %9 . 4f 
>Spv02_f inal, %9 . 4f 
>Spa02_unf i 1 , % 9 . 4 f 
>Spa02_fill,%9.4f 
>Spa02_dual,%9.4f 
DVHS,%9.4f 
DVHS0,%9.4f 
DVAST,%9.4f 
DVAST time,%9.4f 
DWST,%9.4f 
DWST time,%9.4f 
>art_omega_pr , % 9 . 4 f 
>ven_omega_pr , % 9 . 4 f 
>time, %9.4f 

HR Results 

T 3,1,8 
Alg 1 

>hr0 (ir),%9.4f 
>hrl (sig_ref ) , %9.4f 
>hr2 (deriv) , %9.4f 
>hr3 (ir/n),%9.4f 
>hr4 (sig/n) , %9.4f 
>hr5 (spec5) , %9.4f 
>hr6 (ANC),%9.4f 
>hr_final,%9.4f 

Bin Results 



340 



T 2,6,8 
Master Bin 
Bin 1 



Bin 2 
Bin 3 
Bin 4 
Bin 5 

ir_rms,%9.4f 
rd_rras, %9.4f 
>cov,%9.4f 
>arterial sat,%9.2f 
>venous sat , %9 . 2 f 
>seed sat,%9.2f 
>numjpeaks , % 9 . 0 f 
>corr, %9.Sf 

,- Calib Results 

T 4,1,10 
C Results 
ir_trans,%9.4f 
rd_trans, %9.4f 
desired_ir, %9.4f ttiA 
desired_rd, %9.4f mA 
actual_ir,%9.4f mA 
actual_rd, %9 .4f mA 
ir_gain,%9.2f 
rd_gain,%9.2f 
ir_delta,%9.4f 
rd_delta, %9.4f 

display strings 

M 40,0 
Calibrating 
S 52,1 

Masimo Stork 

S 53,0 

Spa02 : 

S 54,0 

Pulse : 

F 14,0 

Tick %12.0f 

F 15,0 

%3.0f %% 

F 1S,0 

%3.0f BPM 

F 17,0 

ADC OOR Count %8.0f 
F 18,0 

Sensor Pos. %8.4f 

Real-Time Plotting 
R 77, 1,400 

Plethysmographies Data (before Masimo) 
Time 

Amplitude 
R 78,1,400 

Plethysmographies Data (after Masimo) 
Time 

Amplitude 
X 12,77,0,9 
IR 

X 13,78,0,12 



end of items 



Exhibit H 



] 



1 /* 

Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 



File : iawa . c 

5 

Description: 

Industry Algorithm Without Averaging 

Public Functions: 
10 init_iawa ( ) ; 

iawa_sat ( ) ; 

Private Functions: 

15 Notes: 

History: 



RJM 11/05/92 

20 RJM 11/20/92 

RJM 12/08/92 

RJM 01/05/93 



Based on s.c (VI. 29) 
use param and result structures 
Make compatible with VSINC 
Compute separate seed for each bin 



# define MODULE_ID 


1015 




#include 


<masimo.h> 


/* 


platform specification 


*/ 


#include 


<memory . h> 


/* 


special memory allocation 


*/ 


# include 


<debug . h> 


/* 


assert, grab 


*/ 


#include 


<math . h> 








# include 


"stork.h" 


/* 


system parameters 


*/ 


# include 


"fir.h" 


/* 


block fir module 


*/ 


# include 


"sategn.h" 


/* 


saturation equations 


*/ 


# include 


"iawa.h" 


/* 


self 


*/ 


# include 


"vector .h" 


/* 


cross corr 


*/ 


#include 


"gddcli.h" 


/* 


General Debug Display Client 




# include 


"spt.h" 









#if (PLATFORM == 1) 

extern VOID 

presentation ( 

FLOAT32 *data, 
INT32 first, 
INT32 num, 
INT32 code) ; 

#endif 
VOID 

init_iawa (void) { 

/* nothing yet */ 



/* this BS is to allow UI . C code to store. 



VOID 

compute_log ( 

FLOAT32 *in, 
FLOAT 3 2 *out, 
INT32 n) { 



/* input data 
/* output data 
/* number of points 



/* This routine computes the natural log of the input signal */ 

/* Because of the filtering in the demodulator, there are times */ 

/* when the resulting signal can go below zero, this log routines */ 

/* limits small values so that the log is gracefully limited. */ 

#define LOG00001 -11.51293 



INT32 



BUG1 (in) ; BUG1 (out) ; BUG1 (n) ; 



for (i = 0; i < n; i++) { 

if (in[i] > 0.00001) out [i] 
else out[i] = LOG00001; 

} 



LOGF(in[i] ) ; 













FLOAT 3 2 


*log_ir, 


/* 


input, normalized ir data 


*/ 


FLOAT32 


*log rd, 


/* 


input, normalized rd data 


*/ 


SAT PARAMS 


*sp. 


/* 


input, parameters 


*/ 


SAT_RESULTS 


*sr) { 


/* 


output , re suit s 


*/ 


INT32 










INT32 


j; 








INT32 


b; 


/* 


bin counter 


*/ 


INT32 


bin_size; 


/* 


point in a bin 


*/ 


INT32 


f_size; 


/* 


size of filtered data 


*/ 


BIN RESULTS 


*r; 


/* 


points to results for each bin 


*/ 


FLOAT32 


*f ir; 


/* 


filtered ir 


*/ 


FLOAT32 


*f_rd; 


/* 


filtered rd 


*/ 


FLOAT32 


*xaxis; 


/* 


for green marker lines 


*/ 



BUGl(log_ir) ; BUG1 ( log_rd) ; BUGl(sp); BUGl(sr); 

assert (log_ir != NULL) ; 
assert (log_rd != NULL) ; 
assert (sp != NULL) ; 
assert (sr != NULL) ; 

f_size = sp->num_points - sp->bp_size; 

bin_size = (f_size - sp->anc_settle_time) / sp->bin_count; 
assert (f_size > 0) ; 

mark_mem() ; 

f_ir = get_mem (f_size) ; 
f_rd = get_mem(f_size) ; 
xaxis = get_mem (1000) ; 

/* band-pass filter IR and RD signals */ 

blkfir (log_ir, sp->num_points, sp->bp_coef, sp->bp_size, f_ir) ,- 

blkf ir (log_rd, sp->num_points, sp->bp_coef, sp->bp_size, f_rd) ; 

/* output analysis data */ 

b=0,- /* build a dataset that demarks the bins divisions */ 

assert (f_size-sp->anc_settle_time < 1000); 
for(i=0; i < f_size - sp->anc_settle_time; i++) 

xaxis [i] =(FLOAT32) (b-5000) ; 
if ( ! (i%bin_size) ) b=10000-b; 



j = sp->anc_settle_time- 1; 
b = f_size - j ; 

GDDCL^CharttSP^AnalysisData.S^, 0.0, 1.0); 
GDDCLI_Data (SPT_AnalysisData,b, &f_rd[j] ) ; 
GDDCLI_Data (SPT_AnalysisData, b, &f_ir [ j ] ) ; 
GDDCLI_Data (SPT_AnalysisData,b, xaxis) ; 

#if (PLATFORM == 1) 

presentation (f_rd, j ,b, 0) ; 
#endif 



calculate correlations 



sr->bins[0] .correlation = VECTOR_CrossCorrO (&f_ir [j ] ,&f_rd[j] ,f_size-j) ; 
for (b = 1; b <= sp->bin_count; b++) 
r = &(sr->bins [b] ) ; 

r->correlation = VECTOR_CrossCorrO (&f_ir [j] , &f_rd[j] ,bin_size) ; 
j+= bin_size; 

/* Square both signals */ 

for (i = 0; i < f_size; i++) 

f_ir[i] *= f_ir[ij; 
f_rd[i] *= f_rd[i] ; 



/* Integrate using binning */ 

sr->bins [0] .ir_rms = 0.0; 
sr->bins [0] .rd_rms = 0.0; 
j = sp->anc_settle_time - 1; 

for (b = 1; b <= sp->bin_count; b++) 

r = &(sr->bins [b] ) ; 

r->ir_rms = 0.0; 
r->rd_rms = 0.0; 

for (i = 0; i < bin size; i++) 
{ 

r->ir_rms += f_ir[j] ; 
r->rd_rms += f_rd[j] ; 
j++; 

} 

/* update master bin */ 

sr->bins [0] . ir_rms += r->ir_rms; 
sr->bins [0] .rd_rms += r->rd_rms; 

r->ir_rms = SQRTF (r- >ir_rms / bin_size) ; 
r->rd_rms = SQRTF ( r- >rd_rms / bin_size) ; 

r->seed_sat = calc_sat (r->rd_rms / r->ir_rms, sp->sat_eqn) ; 



/* compute master bin results */ 
r = &(sr->bins [0] ) ; 

r->ir_rms = SQRTF (r- >ir_rms / f_size) ; /* was bin_size HGK */ 

r->rd_rms = SQRTF (r- >rd_rms / f_size) ; 

r->seed_sat = calc_sat (r->rd_rms / r->ir_rms, sp->sat_eqn) ; 
f ree_mem ( ) ; 



Exhibit I 



X /* 

Copyright (c) Masimo Corporation (tm) 1992. All Rights Reserved. 



File: satfil.c 

Description: Spa02 output filtering algorithms 

Public Functions: 
reset_satf il 
sat_filter 

Private Functions: 
max_dif f_alg 



RJM 
RJM 
RJM 
RJM 
RJM 
RJM 
RJM 



10/09/92 
10/12/92 
11/06/92 
01/16/93 
01/18/93 
01/26/93 
01/26/93 



Initial version 

Allow multiple filtering streams 

Bring up to coding standards 

Merge binning algorithms, DVHS filters, etc 

If COV is low, and 1 peak don't report Sv02 

Limit final Spa02 by highest seed of all bins 

Compare max-diff values to Spa02_final 

*/ 



#define M0DULE_ID 



#include 


<masimo . h> 


/* 


platform description file 


*/ 


#include 


" stork. h" 


/* 


system parameters 


*/ 


# include 


"satfil.h" 


/* 


self 


*/ 


# include 


"sateqn.h" 


/* 


calc_omega 


*/ 


ftinclude 


<memory.h> 


/* 


Special memory routines 


*/ 


# include 


<debug . h> 


/* 


assert , grab 


*/ 


# include 


<vector . h> 


/* 


vector functions 


*/ 


#include 


<math.h> 


/* 




*/ 



VOID 

reset_satf il ( 

SAT_PARAMS *sp, 
SAT_RESULTS *sr) { 

BUGl(sp) ; BUGl(sr) ; 

assert (sp 1= NULL) ; 
assert (sr != NULL) ,- 



sr->DVHS = sp->DVHS_init; 
sr->DVHS0 = sp~>DVHS0_init; 
sr- >Spa02_f inal = 0.0; 
sr->Spa02_unfil = 0.0; 
sr->Spa02_f ill = 0.0; 
sr->Spv02_f inal = 0.0; 
sr- >Spa02_dual = 0.0; 
sr->DVAST = 0.0; 
sr->DWST = 0.0; 

sr->DVAST_time = -sp->DVAST_timer 
sr->DWST_time = -sp->DWST_timer 



/* input, saturation parameters 
/* inout, saturation results 



65 /* 

Name: max_diff_alg 



This routine finds the two bins with the largest difference 
between arterial and venous sat and of these two returns 
the saturation which is closest to the previous sat. For 
more detail see Stork design note 27. 



Select the bins that have a good arterial and t 
saturation and order these bins based on the difference between 
the venous and the arterial saturation. 

Select the two bins with the greatest difference. If there 
are no bins resulting from step 1. 

Of these two bins select the one closest to the previous 
Spa02_final value. 

If there is no previous value, then return the largest value. 



diff alg( 










SAT P ARAMS 


*sp, 


/* 


input, saturation parameters 


*/ 


SAT_RESULTS 


*sr) { 


/* 


inout, saturation results 


*/ 


INT32 


i; 








BIN RESULTS 


*b; 


/* 


general bin ptr 


*/ 


INT32 


*good bins; 


/* 


array of good bin numbers 


*/ 


FLOAT32 


*diffs; 


/* 


array of differences 


*/ 


INT32 


num_good_bins ; 


/* 


number of good bins 


*/ 


FLOAT32 


sat; 


/* 


return value 


*/ 


FLOAT32 


satl; 


/* 


arterial sat of bin with largest 


diff*/ 


FLOAT32 


sat2 ; 


/* 


sat of bin with 2nd largest diff, 


*/ 


BUG1 (sp) ; BUG1 (sr) ; 








assert (sp ! = 


= NULL) ; 








assert(sr != NULL); 









mark_mem ( ) ; 

good_bins = get_mem (sp- >bin_count) ; 
diffs = get_mem(sp->bin_count) ; 

/* select bins that have good arterial and venous sat */ 
/* R154 : bin must also have good correlation HGK */ 

num_good_bins = 0 ; 

for (i = 1; i <= sp->bincount; i++) { 
b = &(sr->bins [i] ) ; 

if (b->arterial_sat > 0 && b->venous_sat > 0 && 
b- correlation >= sp- correlation) 



{ 



good_bins [num_good_bins] = i; 

dif f s [num_good_bins] = b->arterial_sat - b->venous_sat; 
num_good_bins++ ; 



/* Find Two bins with largest difference */ 

switch (num_good_bins) { 

case 0: 

sat = sr->Spa02_unf in- 
break; 



case 1: 



i = max__index (dif f s, num_good_bins) ; 

sat = sr->bins [good_bins [i] ] . arterial_sat ; 

break; 

default : 

/* find bin with largest diff */ 

i = max_index(dif f s, num_good_bins) ; 

satl = sr->bins [good_bins [i] ] .arterial_sat; 

/* find bin with, second largest diff */ 

diffs[i] = 0.0; 

i = max_index(dif f s, num_good_bins ) ; 

sat2 = sr->bins [good_bins [i] ] -arterial_sat; 

/* If there is no previous Spa02 choose highest sat */ 

if (sr->Spa02_f inal ==0.0) { 
sat = satl; 

if (sat2 > satl) sat = sat2 ; 
} else { 

/* Otherwise, choose one closest to Spa02_final */ 
sat = satl; 

if (FABSF (satl - sr->Spa02_final) > 

FABSF (sat 2 - sr->Spa02_f inal) ) sat = sat2; 



f ree_mem ( ) ; 
return ( sat ) ; 

} 



sat_f ilter 

This routine takes the master power curve results and the 
five bins and computes an overall arterial and venous 
saturation based on the following algorithm. For more 
detail see Stork design note 27 . 



#define MAX_SPA02 99 
#define MIN_SPA02 40 

VOID 

sat_f ilter ( 

SAT_PARAMS *sp, 
SATJRESULTS *sr) { 

BIN_RESULTS *mb; 
FLOAT 3 2 delta; 
INT32 i; 
FLOAT 3 2 max_seed; 



/* max sat for final filter output 
/* min sat for final filter output 



/* master bin ptr 

/* maximum seed from all bins 



BUGl(sp) ; BUGl(sr) ; 



assert (sp != NULL); 



assert (sr != NULL) ; 

/* If master bin cov is less than limit then use its sat value */ 
/* otherwise, use the sat of the max_diff_alg of the bins. */ 

mb = &(sr->bins [0] ) ; 

if (mb->cov < sp->cov_limit) { 

sr->Spa02_unfil = mb->arterial_sat; 

/* Update DVHS */ 

if (mb->arterial_sat > sr->DVHS0) { 
sr->DVHS = mb->arterial_sat; 
sr->DVHS0 = sr->DVHS; 

} 

/* DVAST algorithm */ 

if (mb->arterial_sat > 0.0 && mb->arterial_sat < sr->DVHS) { 
sr->DVAST = mb->arterial_sat; 
sr->DVAST time = sr->time; 

} 

} else { 

sr->Spa02_unf il = max_dif f_alg (sp, sr); 
/* DWST algorithm */ 

if (mb->venous_sat > 0.0 && mb- >num_peaks > 1) { 
sr->DWST = mb->venous_sat; 
sr->DWST_time = sr->time; 

, 1 

/* first stage of filtering (DVHS) */ 
sr->Spa02_f ill = sr->Spa02_unfil; 

if (sr->DVHS < sr->Spa02_fill) sr->Spa02_f ill = sr->DVHS; 

/* second stage of filtering */ 
/* If there is motion and a previous value, limit the change */ 
/* to max_sat_delta . */ 

if (mb->cov > sp->cov_limit && sr->Spa02_f inal > 0.0) { 
delta = sr->Spa02_fill - sr->Spa02_f inal; 

if (delta > sp- >max_sat_delta) delta = sp- >max_sat_delta; 
if (delta < -sp->max_sat_delta) delta = -sp->max_sat_delta; 
sr->Spa02_f inal += delta; 
} else { 

sr->Spa02_f inal = sr->Spa02_f ill; 

} 

/* Find largest seed from all bins */ 

max_seed = mb->seed_sat ; 

for (i = 1; i <= sp->bin_count; i++) { 

if (sr->bins [i] .seed sat > max seed) max_seed = sr->bins [i] . seed_sat; 

} 

if (sr->Spa02_f inal < max_seed) sr->Spa02_f inal = maxseed; 
if (sr->Spa02_final > MAX_SPA02) sr- >Spa02_f inal = MAX_SPA02 ; 
if (sr->Spa02_final < MINJ3PA02) sr->Spa02_f inal = MIN_SPA02 ,- 

/* Venous logic */ 

if (mb->cov < sp->cov_limit && mb- >num_peaks < 2) { 
sr- >Spv02_f inal = 0.0; 



} else { 

sr->Spv02_f inal = mb->venous_sat; 

275 } 

/* Compute arterial and venous sat for use by Pulse Rate algorithms */ 

sr->art_omegaj>r = 0.0; 
280 if (sr->time - sr- >DVAST_time < sp->DVAST_timer) { 

if ( FABSF ( sr- >DVAST - sr->Spa02_f ill) < sp->DVAST_delta) { 
sr- >DVAST_time = sr->time; 

} 

sr->art_omegaj>r = calc_omega (sr- >DVAST, sp->sat_eqn) ; 
285 } else { 

if (sr->Spa02_fill > 0.0) { 

sr- >art_omega _pr = calc_omega(sr->Spa02_fill, sp->sat_eqn) ,- 



290 

sr- >ven_omega_jpr = 0.0; 

>time - sr- >DWST_time < sp->DWST_timer) { 
(FABSF (sr - >DWST - mb- >venous_sat) < sp- >DWST_delta) { 
sr->DWST_time = sr->time; 

>ven_omega_pr = calc_omega (sr->DWST, sp- >sat_eqn) ; 

(mb->venous_sat > 0.0 && 

mb->cov > sp->cov_limit && mb- >num_peaks > 1) { 
sr->ven_omega_pr = calc_omega(mb->venous_sat, sp->sat_eqn) ; 



305 

} /* satfil */ 



310 #ifdef STANDALONE 

#include <stdlib.h> 
ftinclude <stdio.h> 

315 VOID main (void) { 



} 

320 #endif 



if (sr- 
295 if 

} 

sr- 
} else 
300 if 



Exhibit K 



1 /* 

Copyright (c) Masimo Corporation (tra) 1992, 1993. All Rights Reserved. 



File: ui.c 



Description: 

This is the STORK software user interface module. This platform 
independant user interface provides the ability to start and stop 
the oximeter, change system parameters, and control the input data. 
There are minor differences in these features listed below: 



Feature 



PC Platform 



DSP Platform 



Start/Stop 
Change system param 
LED calibration 
PAN through data 
Rerun same data 
Standard Mode 

Public Functions: 
init_ui 

user_interface 

Private Functions: 
store_raw 
store_param 
store_results 
single_step 
test_pattern 



Notes : 




History: 




RJM 


12/09/92 


RJM 


12/28/92 


RJM 


12/30/92 


RJM 


12/30/92 


RJM 


01/18/93 


RJM 


01/18/93 


RJM 


01/18/93 


RJM 


01/25/93 


RJM 


01/25/93 


RJM 


01/25/93 


RJM 


02/02/93 


RJM 


02/11/93 


RJM 


02/12/93 


RJM 


02/19/93 


RJM 


02/23/93 


RJM 


02/24/93 


HGK 


03/08/93 



Yes 
Yes 
Yes 



Yes 
Yes 
Yes 



Standard initialization 
Start the user interface 



send raw data via LSL 

send parameters via LSL 

send results via LSL 

Run a single snapshot of data 

send a test pattern via LSL 



Initial Version based on main.c 

New HR module vsi_hr 

New default values for parameters 

Add peak_epsilon parameter 

Add dvhs__init and dvhs0_init parameters 

Add hrfil module 

Keep track of time 

Set default for DVHS0 to 95 

Make SAT egn parameter 

hr_enable param 

Add command for DSP Platform version 
Validate data after get_snapshot 
New default values for coefs 
Integrate file list feature from HK. 
Add Calibration parameters 
Improve ADC Snapshot 
move data_ptr code to usrl_c 



55 



#define M0DULE_ID 


1024 




# include 


<masimo.h> 


/* 


Platform description 


*/ 


#include 


<pckeybrd . h> 


/* 


key definitions 


*/ 


# include 


<memory . h> 


/* 


memory manager 


*/ 


# include 


< debug. h> 


/* 


virtual platform debugger 


*/ 


# include 


<math.h> 








# include 


" sateqn.h" 


/* 


saturation equations 


*/ 


# include 


" stork. h" 


/* 


System parameters 


*/ 


# include 


" adc.h" 


/* 


Demodulation module 


*/ 


iinclude 


"usrl.h" 


/* 


Filtering and Decimation 


*/ 


#include 


"calib.h" 


/* 


LED calibration 


*/ 


# include 


"vsinc .h" 


/* 


VSI noise cancellation module 


*/ 


# include 


"satfil.h" 


/* 


Saturation Filters 


*/ 


#include 


- hrfil. h- 


/* 


Pulse Rate Filters 


*/ 





#include 


"lsl .h" 


/* 


low speed link module 


*/ 


70 


# include 


"iawa.h" 


/* 


Industry algorithm without averaging 


*/ 




#include 


"post.h" 


/* 


post band pass filter 


*/ 




# include 


"vsihr .h" 


/* 


VSI Heart Rate Module 


*/ 




#include 


"vsipl.h" 


/* 


VSINC Clean photoplethysmograph 


*/ 




#include 


" gddcli.h" 








75 


#include 


"spt.h" 


/* 


Stork Packet Types 


*/ 




# include 


"vsidual .h" 


/* 


dual ANC sat version 


*/ 




# include 


"ui.h" 


/* 


Self 


*/ 




#if (PLATFORM == 1) 









#include <stdio.h> 
static BOOL presentation_f lag; 
ttendif 



typedef struct 
FLOAT32 
FLOAT32 
FLOAT32 
FLOAT32 

} PARAMJDESC; 



min_val 
max_val 
inc_val 
def val 



/* minimum value 

/* maximum value 

/* amount to change 

/* default value 



typedef struct { 



SAT P ARAMS 


*vp; 


/* 


pointer to vsinc params 


*/ 


SAT RESULTS 


*vr ; 


/* 


pointer to vsinc results 


*/ 


FLOAT32 




/* 


pointer to power curve 


*/ 


HR PARAMS 


*vhp; 


/* 


pointer to Heart Rate params 


*/ 


HR RESULTS 


*vhr ; 


/* 


pointer to Heart Rate results 


*/ 


FLOAT 3 2 


*f spec; 


/* 


frequency spectrum 


*/ 


FLOAT32 


*pleth; 


/* 


clean plethysmograph 


*/ 


CALIB PARAMS 


*cp; 


/* 


pointer to calibration params 


*/ 


CALIB RESULTS 


*cr; 


/* 


pointer to calibration results 


*/ 



ALGORITHM S PEC ; 



The following constants control the user adjustable 
parameters. 



#define MAX_ALG0R I THMS 
#define NUM_PARAMS 



# define 
#define 
#define 
#define 
# define 
#define 
#define 
#def ine 
#define 
#define 
ftdefine 
#define 
idefine 
# def ine 
# def ine 
tdefine 
#define 
#define 
ttdefine 
#define 
#def ine 
#define 
idefine 
#def ine 
#define 



P_cov_limit 

P_num_point s 

P_anc_type 

P_lambda 

P_min_error 

P_anc_sett le_t im< 

P_num_cells 

P_peak_limitl 

P_peak_epsilon 

P_search_f irst 

P_search_step 

P_search_size 

P_bin_count 

P_sampl ing_rate 

P_dvhs_init 

P_dvhsO_init 

P_DVAST_timer 

P_DVAST_delta 

P_DWST_timer 

P_DWST_delta 

P_max_sat_delta 

P_max_hr_delta 

P_sat_eqn 

P_hr_enable 

P_cp_target 



#define P_cp_error_band 
#define P_cp_settle_time 
#define P_cp_ir_gain 
tdefine P_cp_rd_gain 
#define P_hr_lambda 
#define P_dual_lambda 
#define P_correlation 

static PARAM_DESC params [NUM_PARAMS] 



/* 








mm 




max 




mc 


def */ 


/* 


cov_limit 


0 


*/ 


0 


0, 


200 


0, 


2 


0, 


40.0 


/* 


numjpoints 


1 


*/ 


100 


0, 


1000 


0, 


10 


0, 


650.0 


/* 


anc_type 


2 


*/ 


1 


0, 


2 


0, 


1 


0, 


1.0 


/* 


lambda 


3 


*/ 


o 


0, 


1 


0, 


0 


0001, 


0.75 


/* 




4 


*/ 


0 


0, 


1 


0 , 


0 


0001, 


0.0005 


/* 


anc settle time 5 


*/ 


0 


0 , 


1000 


o , 


1 


o, 


149.0 


/* 


num_cells 


6 


*/ 


1 


o, 


50 


o , 


1 


o, 


12.0 


/* 


peak_limitl 


7 


*/ 


0 


o , 


1 


o , 


0 


01, 


0.02 


/* 


peak_epsilon 


8 


*/ 


0 


o, 


1 


o , 


0 


01, 


0.03 


/* 


search_f irst 


9 


*/ 


0 


o, 


120 


o , 


0 


25, 


104.0 


/* 


search step 




*/ 


0 


o, 


5 


o , 


0 


10, 


0.5 


/* 


search_size 


11 


*/ 


0 


o, 


300 


o , 


5 


o. 


195.0 


/* 


bin_count 


12 


*/ 


1 


0 , 


5 


0, 


- 1 


o. 


5.0 


/* 


sampling rate 


13 


*/ 


10 


0, 


200 


o. 




o, 


62.5 


/* 


dvhs_init 


14 


*/ 


40 


0, 


105 


0, 




1, 


99.0 


/* 


dvhsO init 


15 


*/ 


40 


o. 


105 


o, 


0 


1, 


95.0 


/* 


DVAST timer 


16 


*/ 


0 


o , 


100 


o , 


1 


o, 




/* 


DVAST delta 


17 


*/ 










0 


1, 


2.0 


/* 


DWST_timer 


18 


*/ 


0 


0, 


100 


0, 


1 


0, 


10.0 


/* 


DWST_delta 


19 


*/ 


o 


0, 


100 


o, 


0 


1, 


4.0 


/* 


max_sat_delta 


20 


*/ 


0 


o, 


20 


0, 


0 


1, 


2.0 


/* 


max_hr_delta 


21 


*/ 


0 


o, 


40 


o, 


0 


1, 


4.0 


/* 


sat_egn 


22 


*/ 


0 


o, 


1 


0, 


1 


0, 


1.0 


/* 


hr_enable 


23 


*/ 


0 


o, 


255 


o, 


1 


o. 


251.0 


/* 


cp_target 


24 


*/ 


0 


0, 


1 


0, 


0 


05, 


0.25 


/* 


cp_error_band 


25 


*/ 


o 


o, 


0 


5, 


0 


01, 


0.10 


/* 


cp_settle_time 


26 


*/ 


0 


0, 


300 


0, 


2 


00, 


75.0 


/* 


cp_ir_gain 


27 


*/ 


0 


o, 


4096 


o, 


0 


1, 


2.0 


/* 


cp_rd_gain 


28 


*/ 


o 


o, 


4096 


o, 


0 


1, 


2.0 


/* 


hr_lambda 


29 


*/ 


o 


1, 


1 


0, 


0.0001, 


0.9999 


/* 


dual_lambda 


30 


*/ 


o 


1, 


1 


o. 


0.0001, 


0 .9999 


/* 


correlation 


31 


*/ 


o 


1, 


1 


o, 


0.0001, 


0.9999 



Static FLOAT 3 2 

VOID 

init_ui (void) { 

/* nothing yet */ 

} 



param_val [MAX_ALGORITHMS] [NUM_P ARAMS] ; 



store raw,plethy data to ASCII files to prepare for presentations 



#if (PLATFORM = 



1) 



VOID 

presentation ( 

FLOAT32 *data, 

INT32 first, 

INT32 num, 

INT32 code) 



{ 



static CHAR filename [65] ; 

CHAR str [65] ; 

FILE *fp; 

INT32 i; 



if ( !presentation_f lag) return; 
switch (code) { 

case 0 : printf ( "Enter B/A filename (.B, .A will be added) : "); 

scanf ( "%s" , filename) ,- 

sprintf (str, "%s.B" , filename) ; 

goto part2; 
case 1 : sprintf (str, "%s. A" , filename) ,- 

t2: 

if ( (fp=fopen(str, "w+ n ) )==NULL) break; 
for (i=6;i<num;i++) { 

fprintf (fp, "%f\n»,data [first+i] ) ; 

f close (fp) ; 
break; 



static VOID 

store_raw ( 

SAT_PARAMS 

FLOAT32 

FLOAT32 

INT32 

INT32 



/* input, parameters 



*/ 

/* input, ir datav */ 
*rd, /* input, red data */ 

size, /* input, number of points */ 

time_stamp) { /* input, starting point number */ 



BUGl(sp); BUGl(ir); BUGl(rd); BUGl(size); BUG1 (time_stamp) ; 

#if (PLATFORM == 1) 

mark_mem() ; 

t = get_mem(size) ; 

assert (t != NULL) ; 

j = sp->bp_size/2 + sp->anc_settle_time; /* 1st pt */ 

k = size - sp->bp_size - sp->anc_settle_time; /* #usable pts */ 



+ ) t[i] = 
+ ) t[i] = 



( FL0AT3 2 ) t ime_stamp , 1.0); 



for(i=0;i<j 
for( ,-i<j+k 
for( ,-i<size ; 

GDDCLI_Chart (SPT_UsrlSnapShot, 3, size, 
GDDCLI_Data(SPT_UsrlSnapShot, size, ir) 
GDDCLI_Data (SPT_UsrlSnapShot, size, rd) 
GDDCLI_Data(SPT_UsrlSnapShot, size, t ) 

GDDCLI_Chart (SPT_RawAnalysisData, 2,k, (FLOAT32) time_stamp, 1.0); 
GDDCLI_Data (SPT_RawAnalysisData, k, &ir[j]); 
GDDCL I_Da t a (SPT_RawAnalysisData, k, &rd[j]); 

free_mem() ; 
ttendif 



static VOID 
store_parameters ( 

INT32 num_algorithms) { 



BUG1 (num_algorithms) ; 

for (i = 0; i < num_algorithms ; i++) { 



^ GDDCLI_Data (SPT_SystemParam, NUM_PARAMS , param_val [i] ) ; 

} 

static VOID 
store_results ( 

ALGORITHM_SPEC *alg, 

INT32 num_algorithms) { 

SAT_PARAMS *vp; 

SAT_RESULTS *vr; 

HR_PARAMS *vhp; 

HR_RE SULTS *vhr; 

FLOAT32 *power; 

FLOAT32 *fspec; /* Frequency SPECtrum */ 

FLOAT32 *pleth; /* clean pleth. */ 

INT32 i; 

INT32 j ; 

FLOAT32 temp; 

BUG1 (alg) ; BUG1 (num_algorithms) ; 

mentry("store_results") ; 

/* frequency spectrum */ 

i = 0; 

vp = alg[i] .vp; 
vr = alg[i] .vr; 
power = alg[i] .power; 
vhp = alg[i].vhp; 
vhr = alg[i] .vhr; 
fspec = alg[i] .fspec; 
pleth = alg[i] .pleth; 

/* frequency spectrum */ 

GDDCLI_Chart (SPT_FreqSpec, 1, vhr->num_bins, vhp->min_hr, 

vhr->BPM_per_bin) ; 
GDDCLI_Data (SPT_FreqSpec, vhr->num_b±ns, alg [i] . fspec) ; 

/* power curve */ 

GDDCLI_Chart (SPT_MasterPowerCurve, num_algorithms, vp- >search_size, 

vp->search_f irst, - (vp- >search_step) ) ; 

for (i = 0; i < num_algorithms; i++) { 

GDDCLI_Data ( SPT_MasterPowerCurve , vp->search_size, alg[i] .power); 

for (j =0; j < vp->bin_count; j++) { 

GDDCLI_Chart (SPT_PowerCurvel + j, num_algorithms, vp- >search_size, 

vp->search_f irst, - (vp->search_step) ) ; 

for (i = 0; i < num_algorithms; i++) { 

GDDCLI_Data (SPT_PowerCurvel + j, vp->search_size, 
^ &(alg[i] .power [ (j+1) *vp->search_slze] )) ; 

} 

/* plethysmography signal */ 

GDDCLl_Chart (SPT_CleanPleth, num_algorithms , pleth_size (vp) , 0.0, 1.0); 
for (i = 0; i < num_algorithms; i++) { 

GDDCIiI_Data(SPT_CleanPleth, pleth_size (vp) , alg [i] .pleth) ,- 



#if (PLATFORM == 1) 

presentation (alg [0] .pleth, 0 , pleth_size (vp) ,1) ; 
#endif 



/* store tabular results */ 

for (i = 0; i < num_algorithms; i++) { 

GDDCLI_Data(SPT_SatResults, SAT_RESULTS_SIZE, (FLOAT32 * ) alg [i] . vr) ; 
^ GDDCLI_Data (SPT_HRRe suits, HR_RESULTS_SIZE, (FLOAT32 * ) alg [ij . vhr) ; 

for (i = 0; i <= alg[0] .vp->bin_count; i++) { 

GDDCLI_Data (SPT_BinResults , BIN_RESULTS_SIZE, 
^ ~ ~" (FLOAT32 *) &(alg[0] .vr->bins [i] ) ) ; 

/* master bins' correlation sent separately (times 100) */ 

temp = alg [0] .vr->bins [0] .correlation * 100.0; 
GDDCLI_Data (SPT_correlation, 1, &temp) ; 

mexit("store results"); 

} 

static INT32 
compute_snapshot_size ( 

ALGORITHM_S PEC *alg, 

INT32 num_algorithms) { 

INT32 i; 
INT32 size; 

BUG1 (alg) ; BUG1 (num_algorithms) ; 

/* find largest snapsnot size */ 
size = 0; 

for (i = 0; i < num_algorithms; i++) { 

if (alg[i] . vp->num_j>oints > size) size = alg[i] .vp->num_points; 

return (size) ; 

} 

static VOID 
new_jpatient ( 

ALGORITHM_SPEC *alg, 

INT32 num_algorithms) { 

INT32 i; 

BUG1 (alg) ; BUG1 (num_algorithms) ; 

for (i = 0; i < numalgorithms; i++) { 

reset_satfil(alg[i] .vp, algti] .vr) ; 
^ reset_hrfil(alg[i] .vhp, alg[i].vhr); 



static VOID 
single_step ( 

ALGORITHM_SPEC *alg, /* input, algorithm structure */ 

INT32 num_algorithms, 

INT32 raw_snapshot_size) { 

INT32 i; /* algorithm */ 

FLOAT32 *ir; /* raw ir data */ 





FLOAT32 


*rd; 


/* 


raw rd data 


*/ 


410 


FLOAT32 


* log_ir ; 


/* 


log of ir data 


*/ 




FLOAT32 


*log_rd; 


/* 


log of rd data 


*/ 




FLOAT32 


*adc_snap; 


/* 


snapshot of adc data 


*/ 




SAT_PARAMS 


*vp; 








415 


SAT_RESULTS 












HR PARAMS 


*vhp- 










HR _ RESULTS 


*vhr; 










FLOAT32 


*power; 










FLOAT 3 2 


*f spec; 


/* 


Frequency SPECtrum 


*/ 


420 


FLOAT32 


*pleth; 


/* 


clean pleth. 


*/ 




INT32 


time; 









BUG1 (alg) ; BUG1 (num_algorithms) ; 



mentry ( " single_step " ) ; 
mark_mem ( ) ; 



BUGl(raw_snapshot_size) ; 



/* Allocate data required for current algorithms 



ir = get_mem 

rd = get_mem 

log_ir = get_mem 

log_rd = get_mem 

adc_snap = get_mem 



{ raw_snapshot_si ze ) 
(raw_snapshot_size) 
( raw_snapshot_size ) 
(raw_snapshot_size) 
( ADC_SNAPSHOT_SIZE ) 



ADC_GetSnapshot (adc_snap) ; 
adc_snap[0] = -1.0; 
adc_snap[l] = 1.0; 

GDDCLI_Chart (SPT_ADCSnapShot, 1, ADC_SNAPSHOT_SIZE, 0.0, 1.0); 
GDDCLI_Data (SPT_ADCSnapShot , adc_snapshot_size, adc_snap) ; 

GDDCLI_SendPacket (SPT_ADC_00R, (FLOAT32) ADC_GetOutOf RangeCount () ) ; 

time = get_snapshot (ir, rd, raw_snapshot_size) ; 

/* validate this snapshot of data */ 

for (i = 0; i < raw_snapshot_size; i++) { 

if (ir[i] > 2.0) assert (0) ,- 

if (ir[i] < -2.0) assert(O); 

if (rd[i] > 2.0) assert (0); 

if (rd[i] <-2.0) assert (0); 

} 

store_raw(alg [0] .vp, ir, rd, raw_snapshot_size, time) ; 

compute_log (ir, log_ir, raw_snapshot_size) ,- 
compute_log (rd, log_rd, raw_snapshot_size) ; 

for (i = 0; i < num_algorithms ; i++) { 

/* setup pointers for current algorithm */ 

vp = alg[i] .vp; 
vr = alg [i] .vr; 
vhp = alg[i] .vhp; 
vhr = algti] .vhr; 

/* allocate variable sized results arrays */ 

alg [i] .power = get_mem (vp->search_size * (vp->bin_count + 1) ) ; 
alg [i].f spec = get_mem(vsi_spectral_size (vhp) ) ; 
alg [i] .pleth = getmem (pleth_size (vp) ) ; 

power = alg [i] .power; 



fspec = alg [i] . fspec; 
pleth = alg[i] .pleth; 

/* time stamp data */ 

vr->time = (FLOAT32) time / vhp->raw_sampling_rate,- 

iawa_sat (log_ir, log_rd, vp, vr) ; 

vsinc(ir, log_ir, log_rd, vp, vr, power); 
sat_f ilter (vp, vr) ,- 

vsi_hr( log_ir, log_rd, vp, vr, vhp, vhr, fspec); */ 

vsi_hr(ir,log_ir, log_rd, vp, vr, vhp, vhr, fspec); 

hr_f ilter (vhp, vhr, vp, vr); 

vsincjpKir, log_ir, log_rd, vr->art_omega_pr , vp, pleth); */ 
vsinc_pl(ir, log_ir, log_rd, 

vr- >art_omega_pr , 
vhp->hr_lambda, /* HR's lambda for plethy. */ 

vp, pleth) ,- 

vr->Spa02_dual = 0.0; 

VSIDUAL_sat (vr- >art_omega jpr , 
vp- >dual_lambda, 
vp, 

ir, rd, log_ir, log_rd, 
&vr->Spa02_dual) ; 

GDDCLI_SendPacket (SPTtimetick, (FLOAT32) time) ; 
GDDCLI_SendPacket(SPT_sat_final, vr->Spa02_f inal) ; 
GDDCLI_SendPacket (SPT_hr_f inal , vhr- >hr_f inal ) ; 

} /* for */ 

store_j?arameters (num_algorithms) ; 
store_results (alg, num_algorithms) ; 

GDDCLI_NewLine ( ) ; 

GDDCLI_FlushBuf fer () ; /* send data */ 

f ree_mem () ; 
mexit("single_step") ; 



static VOID 
test_pattern (void) { 

/* send a test pattern to verify graphing stuff */ 

} 

static VOID 
set_alg ( 

ALGORITHM_SPEC *alg, 

INT32 num_alg) { 

SAT_PARAMS *vp; 

HR_PARAMS *vhp; 

CALIB_PARAMS *cp ; 



BUG1 (alg) ; BUG1 (num_alg) ; 
for (i = 0; i < num_alg; i++) 



vp = alg [i] . vp; 



vp->cov_limit 

vp- >numjpoints 

vp->anc_type 

vp->lambda 

vp- >min_error 

vp- >anc_settle_time 

vp- >num_cells 

vp- >peak_l imitl 

vp- >peak_epsilon 

vp->bp_coef 

vp->bp_size 

vp- >search_f irst 

vp- >search_step 

vp- >search_size 

vp- >bin_count 

vp->DVHS_init 

vp- >DVHSO_init 

vp- >DVAST_timer 

vp- >DVAST_delta 

vp- >DWST_timer 

vp- >DWST_delta 

vp- >max_sat_delta 

vp->sat_eqn 

vhp = alg[i].vhp; 

vhp- > resolution 

vhp- >raw_sampling_rate 

vhp->min_hr 

vhp->max_hr 

vhp->nst 

vhp- >max_hr_delta 
vhp- >hr_enable 

vhp- >hr_lambda 
vp- >dual_lambda 
vp- ^-correlation 

cp = algfi] .cp; 
cp->rd_min_mA 
cp- >rd_max_raA 
cp - > i r_in injmA 
cp->ir_max_mA 
cp->target 
cp->error_band 
cp - > i r_dr ive_l 
cp- >rd_drive_l 
cp- >ir_drive_2 
cp- >rd_drive_2 
cp- >settle_time 
cp->ir_gain 
cp->rd_gain 
cp->mu 



= param_val[i] 
= (INT32)param 
= ( INT3 2 ) param 
= param_val [i] 
= param_val [i] 
= (INT32) param 
= (INT32)param 
= parani_val [i] 
= param_val [i] 
= post_coef ; 
= post_size; 
= param_val [i] 
= param_val [i] 
= (INT32)param 
= (INT32) param 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val ti] 
= param_val [i] 
= (INT32)param 



[P_cov_limit] ; 
i_val [i] [P_num_points] ; 
_val [i] [P_anc_type] ; 
[P_lambda3 ; 
[P_min_error] ; 

_val [i] [P_anc_settle_time] ; 
_val [i] [P_num_cells] ; 
[P_peak_limitl] ; 
[P_jpeak_epsilon] ,- 



[P_search_f irst] ; 
[P_search_step] ; 
_val [i] [P_search_size] ; 
_val [i] [P_bin_count] ; 
[P_dvhs_init] ; 
[P_dvhsO_init] ; 
[P_DVAST_timer] ; 
[P_DVAST_delta] ; 
[P_DWST_timer] ; 
[P_DWST_delta] ; 
[P_max_sat_delta] ; 
i_val [i] [P_sat_eqn] ; 



= 1.0; 

■■ (INT32)param_val [i] [P_sampling_rate] ; 

■ 25.0; /* bpm */ 
■■ 225.0; /* bpm */ 

■■ 0.04; /* thresshold */ 

■ param_val [i] [P_max_hr_delta] ; 

: (lNT32)param_val [i] [P_hr_enable] ; 

= param_val[i] [P_hr_lambda] ,- 

■ param_val [i] [P_dual_lambda] ; 

■ param_val [i] [P_corr elation] ; 



20.0; /* mA */ 

50.0; /* mA */ 

5.0; /* mA */ 

50.0; /* mA */ 

param_val [i] [P_cp_target] ; 
param_val [i] [P_cp_error_band] ; 

/* mA */ 
/* mA */ 
/* mA */ 



5.0 
> 20.0 
10.0 
: 25.0 



/* mA 



= (INT32) param_val [i] [P_cp_settle_time] ; 
= param_val [i] [P_cp_ir_gain] ; 
= paramjval ti] [P_cp_rd_gain] ; 

: 0.1; 



#define DATA_STEP 57 /* amount to scroll raw data */ 
VOID 

user_interf ace ( 

INT32 num_algorithms) { 

ALGORITHMS PEC *alg; 

INT32 alg_num; /* current algorithm 





INT32 


param_num; 


/* 


current parameter number 


*/ 


615 






/* 


quit program 


*/ 




BOOL 


automation; 


/* 


quit automated session? 


*/ 




INT32 


lsl_cmd; 


/* 


Low Speed Link command 


*/ 




FLOAT32 


*ir; 


/* 


for data preview 


*/ 


620 


FLOAT3 2 


*rd; 










INT32 


preview size = 


650 








INT32 


ss_size; 


/* 


snapshot size 


*/ 


625 


INT32 


i; 










INT32 


j; 










BOOL 


cal ok; 










INT32 


time ; 










FLOAT32 


f ; 









630 

BUG1 (num_algorithms) ; 



assert (num_algorithms <= MAX_ALGORITHMS) ; 

635 mark_mem(); 

alg = get_mem(num_algorithms * sizeof (ALGORITHM_SPEC) ) ; 

/* now allocate structures for each algorithm */ 

640 for (i = 0; i < num_algorithms ; i++) { 

alg[i].vp = get_mem( sizeof (SAT_PARAMS) ) ; 

alg[i] .vr = get_mem (sizeof ( SAT_RESULTS ) ) ; 

alg[i] .vhp = get^mem (sizeof (HRJ?ARAMS) ) ; 

alg[i].vhr = get_mem (sizeof (HR_RESULTS) ) ; 
645 alg[i].cp = get_mem (sizeof (CALIBJ? ARAMS) ) ; 

alg[i].cr = get_mem (sizeof ( CALIBRESULTS ) ) ; 



/* init param structures */ 

650 

for (i = 0; i < num_algorithms; i++) { 
for (j = 0; j < NUM_PARAMS ; j++) { 

param_val [i] [j] = params [ j ] .def_val ; 

655 } 

param_num = 0; 
alg_num = 0; 
quit = FALSE; 

660 

ir = get_mem (preview_size) ; 
rd = get_mem(preview_size) ; 

new_file (0) ; 

665 

setalg(alg, num_algorithms) ; 
new_patient (alg, num_algorithms) ; 

store_parameters (num_algorithms) ; /* display initial parameters */ 
670 GDDCLI_FlushBuf f er ( ) ; 

CALIB_Init(alg[0] .cp, alg[0].cr); 
do { 

GDDCLI_SendPacket (SPT_CalibrateMsg , 1.0); 
675 GDDCLI_FlushBuffer () ; 

cal_ok = CALIB_TwoPointCal(alg[0] .cp, alg[0].cr); 
GDDCLI_Data (SPT_CalibRe suits , CALIB_RESULTS_SIZE, 

(FLOAT 3 2 *) (alg[0] .cr) ) 

} while (cal_ok == FALSE) ; 
680 GDDCLI_SendPacket(SPT_CalibrateMsg, 0.0); 



GDDCLI_FlushBuffer() ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp, ir, rd, preview_size, time); 
GDDCLI_FlushBuffer() ; 

presentation_f lag = FALSE; 

while (quit == FALSE) { 

GDDCLI_HighLight(SPT_SystemParam, param_nutn, alg_num, 

param_val [alg_num] [param_num] ) ; 

GDDCLI_FlushBuffer() ; 
lsl_cmd = lsl_read(); 

GDDCLl_Normal (SPT_SystemParam, param_num, alg_num, 

param_val [alg_num] [param_num] ) ; 

GDDCLI_FlushBuffer<) ; 

switch (lsl_cmd) { 

case NLT: /* left arrow */ 
if (alg_num > 0) { 
alg num -= 1; 

} 

break; 

case NRT: /* right arrow */ 

if (alg_num < num_algorithms-l) { 
alg_num += 1; 



case NUP: /* up arrow */ 
if (param_num > 0) { 
param_num -= 1; 



case NDN: /* down arrow */ 

if (param_num < NUM_PARAMS-1) { 
param_num += 1; 



case NINS : 

725 f = params [param_num] . inc_val * 100.0; /* HGK */ 

goto kplus ; 

case NPLUS: 

f = params [param_num] . inc_val ; /* HGK */ 

73 0 #if (PLATFORM == 1) 

if (shift_key() ) f*= 10.0; 

#endif 
kplus : 

735 param_val [alg_num] [param_num] += f; 

if (paramval [alg_num] [param_num] > params [param_num] .max_val) 
param_val [alg_num] [param_num] = params [param_num] . max_val ; 
break ; 

740 case NDEL : 

f = params [param_num] . inc_val * 100.0; /* HGK */ 

goto kminus; 

case NMINUS: 

745 f = params [param_num] . inc_val ; /* HGK */ 

#if (PLATFORM ==1) 

if (shift keyO ) f*= 10.0; 

ttendif 



param_val [alg_num] [param_num] -= f; 

if (param_val [alg_num] [parara_num] < params [paramjnum] .min_yal) 
param_val [alg_num] [param_num] = params [param_num] . min_val ; 
break ; 



case NRTUP: /* page up */ 

move_snapshot ( -DATA_STEP) ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp, ir, rd, preview_size, time); 

GDDCLI_FlushBuf f er ( ) ; 

break; 

case NRTDN: /* page down */ 

move_snapshot ( DATA_STEP) ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp, ir, rd, preview_size, time) ; 

GDDCLI_FlushBuffer () ; 

break; 

case SCR: /* continuous operation */ 
set_alg(alg, num_algorithms) ; 

ss_size = compute_snapshot_size(alg, num_algorithms) ; 
* calibrate2 (system_params.cal_target) ; */ 
reset_file_num() ; 
new_file(0) ; 

automation = TRUE; /* 'continue automated session' */ 

do 

{ 

time=0; 

new_patient (alg, num_algorithms) ; 

do 

{ 

single_step (alg, num_algorithms , ss_size) ; 
move_snapshot (USR1_AUT0_STEP) ; 

/* mating get_snapshot ( ) is in single_step . . */ 
if (lsl_status() ) 

automation = FALSE; /* discontinue session */ 

} 

while (automation && !end_of_data ( ) ) ,- 
if (automation) 

GDDCLI NewFileO ; 

} 

while (automation && ! end_of_f iles ( ) ) ; 
break; 

case SCS: /* single step */ 

set_alg(alg, num_algorithms) ; 

ss_size = compute_snapshot_size (alg, num_algorithms) ; 

single_step (alg, num_algorithms, ss_size) ; 

break; 

"C calibrate command */ 

case SCT: /* "T" test pattern command */ 
test_jpattern( ) ; 
break; 

case NESC: /* escape */ 
quit = TRUE; 
break; 

case SCF: /* "F" new file command */ 
new_file(l) ; 



time = get_snapshot (ir, rd, preview_size) ; 
store_raw(alg[0] .vp,ir, rd, prev±ew_size, time) ; 
GDDCLI_FlUShBUf f er ( ) ; ; 

break; 

case SCB: /* "B" back file */ 
new_f ile (-1) ; 

time = get_snapshot (ir, rd, preview_size) ; 
store_raw(alg [0] .vp, ir, rd, preview_size, time); 
GDDCLI_FlushBuf f er ( ) ; 
break ; 

case SCH .- 

param_val [algjnum] [P_num_cells ]=2; 
param_val [alg_num] [P_search_size] =25; 
break; 

#if { PLATFORM == 1) 

case SCW : presentation_f lag = 1- presentation_f lag; 
break; 

#endif 

default: /* ignore unknown commands */ 
break; 

} /* switch */ 

} /* while */ 

f ree_mem ( ) ; 
} /* user_interface */ 



